Ion-acoustic solitons in warm magnetoplasmas with super-thermal electrons 
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In this work, the phenomenon of formation of localised electrostatic waves (ESW) or soliton is 
considered in a warm magnetoplasma with the possibility of non-thermal electron distribution. The 
parameter regime considered here is relevant in case of magnetospheric plasmas. We show that devia- 
tion from a usual relaxed Maxwellian distribution of the electron population has a significant bearing 
in the allowed parameter regime, where these ESWs can be found. We further consider the presence 
of more than one electron temperature, which is inspired by recent space-based observations^. 

I. INTRODUCTION 

Electrostatic solitary waves (ESW) or solitons are common occurrences in the near-earth plasmas and are routinely 
observed in the boundary layers and turbulent regions of the magnetosphere^. In recent years, many authors have 
tried to model these ESWs with the help of different physical models^—. Most of the observational and theoretical 
studies have indicated that these ESWs are basically potential structures (compressive and rarefactive) with density 
structures (enhancement and compression) and weak double layers 1 .. Solitons or double layers are plasma sheaths 
(discontinuities) moving in a plasma, which require competing species of charged particles with different masses 
(inertia) and charges. In an electron-ion plasma, these requirements are fulfilled due to colder ions with their large 
inertia, which helps in formation of plasma sheaths. So, an electron sheath in a plasma can be effectively modified by 
equally hot ions 8 -. 

However, majority of these theoretical analyses do not consider the effect of magnetic field into account. A full- 
blown analysis of solitary waves where magnetic perturbation is self-consistently considered can be very rigorous 9 . 
Several space-borne experiments have observed ESWs moving parallel to the ambient magnetic field in various near- 
earth plasma environments such as the solar wind, magnetosheath and magnetotail regions, and auroral zone—— . 
These ESWs are, in general, bipolar structures moving in the direction of the background magnetic field 1 —. In the 
auroral region, these ESWs are reported as negative potential structures travelling upward along the auroral magnetic 
field line o 15 ' 16 . These observations are also supported by those from the Freja satellite dat a 17 ' 18 . Besides, space 
plasmas, in general, can be largely modeled with Maxwellian velocity distribution. However, advancement of satellite 
based technologies in recent years, has led to the realization that most of these plasmas, especially the near-earth 
plasmas, have high energy tails and heat-flux shoulders, which may be attributed to the fact that these plasmas are 
quite inhomogeneous and semi-collisionlessi 9 - - — . Subsequently, it has been established that these plasmas are best 
modeled by a generalized Lorentzian or kappa distribution (especially the electron distributions) rather than by a 
pure Maxwellia n 22 ' 23 . Experimental observations on solar wind plasmas, in recent years, have established that these 
plasmas have a spectral index 24 ~ 2.8. In the Earth's magnetosphere, the spectral index is typically observed 2 ^ in the 
range 2 < k < 8. We, in this Chapter, consider these ESWs in the presence of a super thermal electron component. 

We, in this work have considered formation of ion-acoustic electrostatic solitary waves in the presence of a back- 
ground magnetic field. We also consider the electron polulation to be super-thermal, modeled through a Lorentzian 
or kappa velocity distribution. We further consider the effect of two species of electrons with different temperatures. 
Presence of multi-temperature electrons in magnetosphere is reported experimentally 10 . In Section I, we present the 
plasma model that we have considered for formation of these ESWs in a warm magnetoplasma. In Section II, we 
consider the solitary wave solutions of the model, where we have incorporated the effect of super-thermal electrons. 
In Section III, we consider two independent components of super-thermal electrons with two different temperatures. 
We however have shown that the two-temperature electrons have only marginal effect on the structure and param- 
eter regimes of the ESWs. Finally, we summarise our conclusions in Section IV. The parameter regime considered 
in this work and the related results can be relevant in explaining large-amplitude ESWs observed in the earth's 
magnetosphere. 

II. BASIC MODEL OF PLASMA 

Below we write down the basic governing equations for a thermal plasma, which is immersed in an external magnetic 
field B , 
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(1) 

— Vp+(wxf2), (2) 

rrii 

(3) 
(4) 

where Eqs.([T]) and represent the continuity equation and conservation of momentum for the ions. The function 
F(<j)) represents equilibrium electron (and ion) density according to the particular velocity distribution function (VDF) 
i.e. Boltzmanian or kappa distribution. The last equation is the equation of state. In the above equations, n and 
n e are the ion and electron densities and f2 = eBo/(rriic) is the ion gyro-frequency. Other symbols have their usual 
meanings and quasi-neutrality is assumed all throughout. We assume the external magnetic field to be in the z 
direction. 

We assume an arbitrary electrostatic perturbation in time and space and define a co-moving coordinate rj — 
l x x + l z z — Vhft, where l x<z are direction cosines and thus defined by the relation l x + l z = 1, and vm is the velocity of 
the nonlinear wave. Far away from the perturbation we assume everything to be stationary and define the boundary 
condition as at r\ — > oo, n —> uq, <\> — > 0, and v — > 0. Without any loss of generality, we can assume that the physical 
quantities are constant along y direction. 



— +V- (nv) = 0, 
at 
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n = n e = F (</)), 

p cc (nrrii) 1 , 



A. Reduction of equations 

We now describe a general procedure for reducing Eqs.([T]|3]). From the continuity equation, we can write 

l x v x + l z v z -v M = —vm — ■ (5) 

n 

From the x, y, and z components of the momentum equation, we can write, 

- v M —v' x = -l x [f{ri) + g(n)] + v y Q, (6) 
n 

-vm — v = -vjl, (7) 
n " 

-v M -v' z = -lAf{n) + g{n)l (8) 
n 

where the (') refers to derivative with respect to the scaled coordinate rj, and f(n) and g(n) functions of n, 

T 

•^( n ) = 1 2 M/3 n » ( 9 ) 

g{n) = —</>', (10) 

rrii 

Here, we have expressed the equilibrium presurre "fPo/n — T, with 7 being the ratio of the specific heats taken as 
5/3. Note that all throughout the calculations, temperature T is expressed in energy units. We now take a derivative 
of Eqs.© and ((HJ with respect to r\ to get the following relations, 

-vm^U--v' x ) = -l x [f'(n) + g'{n)] + v x - — , (11) 
n \ n J n v M 

-vm- f< - -v' z ) = -h[f'(n)+g'(n)}, (12) 
n \ n J 

where we have substituted for v' y from Eq.®. By differentiating Eq.©, successively with respect to 77, we get, 

l x v' x + l z v' z = h(n) = v M ^n', (13) 

1 II , 1 II 1 \ o n 1 l\2 n II / -1 t\ 

lxV x + l z v z = p(n) = -2v M ^(n ) - v M —^n . (14) 
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Using Eas. (|imi4|) . we get, 
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where we have used the condition l 2 + l 2 z = 1 and substituted for l x v x from Eq.([5]). Eq.© can be integrated to have, 



v M n 



n[f(n) + g{n)} dn + C\ = q(n), 



(16) 



where C\ is the constant of integration, to be evaluated by imposing the boundary conditions. So, finally, using 
Eq. (ll6l) . from Eq. (ll5|) . we arrive at a single nonlinear, second order differential equation for the system as, 

-lf{n)+g'{n)]+tf f " s ' " 
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l 7 — —q(n). 

n v M 



(17) 



We note that the arbitrary functions /, g, h,p, and q can be written entirely as functions of n. 



III. SOLITON SOLUTIONS 



In this section, we consider, in general, the electrons to be super thermal governed by Lorentzian or kappa velocity 
distribution function (VDF). The density of super thermal electrons with a kappa VDF can be written as^r— , 

1/2 -re 



2c c 



n = no yl — 

which reduces to the familiar Boltzmann relation in the limit k — > 00. The temperature is represented by 
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We note that validity of this particular VDF requires that the spectral index k > 3/2. The expression for density 
[Eq. (fT8|) ] can be inverted to express the potential <j> as, 



1 TOp 



1- — 



(20) 



Using the above expression for <p in the expression for gin), and carrying out the integration in Eq. (|16p . we finally 
write down the explicit nonlinear differential equation in n from Eq. (|17p as, 



where 
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(21) 

(22) 

(23) 
(24) 



The above equations i.e. Eqs. (|2H24[) are written in terms of normalized variables. The ion density n is normalized to 
its equilibrium density no, velocities to the ion-sound velocity v s — \JT\Jrn~i, potential 4> to T e /e, and length to the 
ratio v s /Vt. The ratio of the ion temperature to the electron temperature is denoted by a — T/T e . Note that with 
this normalizations, the normalized ion velocity Vm is the Mach number. 

Equation (|2ip describes the behavior of nonlinear ion-acoustic wave including solutions of solitary waves. However, 
in order to analytically isolate existence of solitary wave solutions, we need to re-cast Eq. ([2~T|) in a form, 

2 



1 f dn 



5U 1 +V(n)= °' 



(25) 



where V(n) is the pseudo or Sagdeev potential 2 —. A potential-well structure of V(n) ensures the existence of solutions 
of solitary waves. The Eq. (|2ip , however, in its present form with arbitrary k can not be reduced to a form represented 
by Eq. (|25p . except in two limiting cases k 00 and k — > 3/2. The first is the well known case of Maxwellian electrons 
and the second is the extreme limit of super thermal electrons. In what follows, we try to use this two limiting cases. 
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Figure 1: (a) Sagdeev potential for Maxwellian electrons and (b) for highly super thermal electrons (k ~ 3/2). The parameters 
are vm = 0.78, l z = 0.7, and a = 0.2. Note the increase in depth and width of the potential in case of super thermal electrons. 

A. Limiting cases 

1. Maxwellian electrons (k — > 00) 

This is a well known case in which Eq. (|21l) can be reduced to the standard form of Eq. (l2~5j) . with the quasi potential 
K(n)as21, 

Voo(n) = [50v 2 M (an 8 ^ + n 2 - u^)]" 1 x (n 6 l 2 z [5(l - n) + 3o"(l - n 5 / 3 )} 2 

~ 5n 5 ^, f [{10(l - n){l 2 z - n) + ncr{9 + 6n 5/3 - 15n 2/3 ) 
+ l 2 z a(6 + 9n 5/3 - 15n)} 

+ 10{l 2 z - l)nlnn] + 25n 4 v 4 M [n - l) 2 )) . (26) 
In Fig[TJa), Sagdeev potential for Maxwellian electrons is shown for a certain range of parameters. 

2. Highly super thermal electrons (k — > 3/2) 

We now carry out an asymptotic expansion of Eq. (|2ip in the neighborhood of n ~ 3/2. After some simple and 
straight forward but lengthy algebra, we can reduce Eq. ([2~T1) to the following form, 



— G{n) = H(n), 



where, 



G{n) = v 2 M \ + (3-2n)- + 3o-n 2/3 , 



H{n) = 2(n- 1)+ — 
vm 



(3 - 2k)ti Inn - -n(n 5/3 - 1) 
5 



Multiplying both sides of Eq. ([27|) by dG(n)/drj and integrating, we to get, 



(27) 

(28) 
(29) 

(30) 



from which, we can finally write the standard form for the quasi potential equation in the form of Eq. (|25l) with the 
quasi or Sagdeev potential as, 

V 3/2 (n) = [50v 2 M {2v 2 M + n(3 -2k- 2n 5 / 3 cr)}]~ :L x (n 6 l 2 z [6a(n 5 ^ 3 - 1) 
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Figure 2: (a) Domains of existence of solitons in the two limiting cases in the (vm, lz,c) parameter space. Solitons exist only 
left sides of the surfaces. The surface on the front (red) is for k ~ 3/2 and the other (blue) is for k —¥ oo. (b) Projection of 
the surfaces showing the regions of existence of solitons in the (vm,I z ) space. The right side boundaries of the regions are for 
higher a. 

- 5 lnn(3 - 2k)] 2 + lOOn 4 ^,/ (n - l) 2 

+ n 4 v|f[100n(3- 2n){(l 2 z - l)(n- 1) - (l\ - n) Inn} 

- 60na{3n + n 5/3 (2n - 5) + l 2 z (2 - 5n + 3n 5/3 )}]) . (31) 

The constant of integration C<x is evaluated by demanding the boundary condition V( n )ln=i = 0. The Sagdeev 
potential for this limiting case is shown in Figfjjb). In Fig Ufa), we have shown the domain of existence of a soliton 
in the (vm — lz — f) parameter space for both these limiting cases. A soliton, if any, can exist only left of the surfaces 
(please see the captions for details). The projection of the domains in the (vm — 1%) space is shown in FigJ^b). The 
left side of the regions are limited by a —¥ 0. In the overlapping region, solitons in both extreme cases can form. 

3. Limiting Mach numbers 

The limiting Mach numbers can be found out by imposing the condition that V(n) must have a local maximum at 
n = 1, which translates to the conditions 

Wl + cr < v M < Vl + <7, (32) 

for Maxwellian electrons and 

l z a/k- 3/2 + ct < v M v/k-3/2 + ct. (33) 

for highly super thermal electrons (k ~ 3/2). The latter condition reduces to l z y[5 < Vm < \fo for k ^ 3/2 . So, 
with highly Lorentzian electrons, solitons can form even with lower values of Mach numbers. In FigOl we have shown 
the dependence of soliton amplitude A and width A, respectively, with the Mach number for both these cases. Note 
that the amplitude A is determined by first zero of the Sagdeev potential away from n = 1 and the soliton width 
A = A/y/d, where d is the maximum depth of the Sagdeev potential determined by the condition dV (n) / dn\ n ^i = 0. 
As can be seen from FigJ31 in case of highly super-thermal electrons, the soliton can be very steep (smaller width) 
and can reach high amplitude in comparison with thermal electrons. 

Note that the range of the valid Mach number is exclusively dictated by the direction of propagation of the soliton 
with the higher limit parallel to the magnetic field (l z — > 1). 

B. Arbitrary k — numerical solution 

As mentioned earlier, Eq. (|2ip can not be reduced to standard form for arbitrary k, which can analytically demon- 
strate the existence of a suitable quasi potential and thus solitary wave solutions, which however, does not rule out 
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Figure 3: (a) Dependences of soliton amplitude A on Mach number, (b) Dependences of soliton width A on Mach number. 
The arrows show the position of lower limits of the respective Mach numbers. The blue (dashed) curves are for k ~ 3/2 and 
the solid (black) curves are for k — > oo (Maxwellian electrons). The parameters are same as in FigfTJ 



existence of soliton-like solutions for arbitrary k. In order to solve the second order, nonlinear differential equa- 
tion, Eq. (|21|) numerically for solitary wave solutions, we use the shooting method, imposing the boundary conditions 
n(i])\r)^-±M = 1 and start with an initial condition 77/(77) |, ;= o — 0, where M is a large number. In practice, we need 
to solve only one half on the 77-axis in the range 77 £ [0, +M] as the the solutions 71(77) are always symmetric in the 
range [-M, +M], which can be easily seen from the invariance of Eq. (f2"Tj) for 77 — > —7;. In FigHJa), we have shown the 
evolution of the soliton as obtained from numerical solution of Eq. (|2ip for arbitrary k. The corresponding bi-polar 
electric field representing the solitary density structure is shown in FigfJJb). Note that the corresponding electric field 
for a density soliton can be found from the relation, 

£?(„) = -I ^. (34) 
n drj 

We note that in the auroral region, the ion to electron temperature ratio er ~ 0.25 and the ion gyro frequency 
(0/27r) ~ 100 HaS 8 .. Consider now various bi-polar electrostatic structures observed in different regions viz. auroral 
region of the ionosphere and earth's magnetosphere. These structures show peak-to-peak variation of electric field 
ranging from few hundreds of micro volts per meter to milli volts per meter within an interval of micro seconds to 
milli seconds. Usually the amplitudes of pulses in the auroral region are larger (milli volts per meter) having a larger 
duration of the order of milli seconds- 1 ^-. Those observed in the magnetosheath region are of much smaller amplitudes 
(micro volts per meter) with a very short duration (micro seconds )*. Note that the bi-polar structures represented 
in FigQJh) are in the rest frame of the solitons, which can be transformed to the duration of the pulses in the rest 
frame of the detector (on board the satellite). In Fig|5J we show various bi-polar electrostatic structures for these 
parameters in real-time units. In all panels of FiglSJ the pulse is of shorter duration and larger height for highly 
Lorentzian electrons than Maxwellian electrons. The peak-to-peak electric field variation due to Lorentzian electrons 
in panel (d) is ~ 250mV/m with a duration ~ 20 ms, which is typical in these region n 15 ' 28 . For the same parameters, 
the peak-to-peak variation is much smaller ~ 90mV/m with a longer duration ~ 40 ms for Maxwellian electrons. The 
pulses due to Maxwellian electrons are comparable to Lorentzian electrons in both pulse height and width only when 
the pulses become very large, which are however, not observed in these regions of space plasma. So, we conclude 
that the deviation of electron temperature from thermal behavior is an essential factor, which needs to be taken into 
account in order to explain the relatively small-amplitude and steep (smaller width) bi-polar electrostatic structures 
in the auroral regions. 

We further note that the effect due to super-thermal electrons is more pronounced when l z ~ 1 or 6 ~ 0°) where 9 
defines the angle between the ambient magnetic field and the direction of propation of the soliton. This is consistent 
with observvation of these solitons in magnetospheric plasmas^*— . 
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IV. TWO-TEMPERATURE ELECTRONS : SMALL AMPLITUDE SOLITONS 



In this section, we consider the physical scenario with two-temperature electrons, which is motivated by the exper- 
imental observations such as those by the Cluster spacecrafts^. We consider two species of electrons — hot and cold 
one. We still assume that both species of electrons are Lorentzian, which gives us freedom to investigate the situation 
both for thermal and super-thermal cases. All mathematical analyses remaining same, differing only in the expression 
denoting the combined density of the electrons [see Eq. lfTS)) ]. 



n e =n c + n h — n c0 ^1 
and the quasi- neutrality condition [see Eq.([3])], 



2e ( 



■2 I + 7l h0 [ 1 



rih, 



2e<t> 
im e B\ 



1/2-K 



(35) 



(36) 



where the subscripts 'c' and 'ft' refer to the cold and hot electrons with 8 c h respectively being functions of T c /,, 
the temperatures of cold and hot electrons. The primary drawback in this formalism that we can not find an 
equivalent Sagdeev potential for solitons of arbitrary amplitudes as Eq. (|35|) can not be inverted to find a unique 
analytical expression for the potential <fi in terms of the densities. In what follows, we shall expand Eq. (|35j) around 
= 0, assuming the analysis for only small-amplitude solitons, so that we can find out the Sagdeev potential which, 
however, will be valid only in the limit of small amplitude. 

A Taylor expansion of Eq. (|35| around <f> = 0, to the zeroth order enables us to write the plasma potential (un- 
normalized) as, 



T c T h ( n 



n c0 - n h0 



\n c0 T h + n h0 T c J \ 2k - 1 



2k- 3 



(37) 



As before, we normalize the variables and write the normalized electron densities as v c ^ — n C Q t f l o/ n o with the condition 
v c + Vh = 1. The temperature ratio of the cold and hot electrons is j3 = T c /Th and a = T/T^. Proceeding as before, we 
can reduce the equations to the form given by Eq. (|25p and write the Sagdeev potential (for small amplitude, n — > 1) 
as 



Vn 



* 4 [100 ffvffSn 2 



20nf lV 2 M {5f 2 (n + l 2 z )Sn 2 + 5n 2 }} 



200^/1^/3 - hv 2 M + f 2 n 3 ) 



(38) 



where Sn 



1 and other quantities are given by, 
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= (2k -1) {l-v h (l ~ 


/?)}, 


h 


= (2k -3)/3, 




8m 


= 5f 2 lz(n 2 - 1) + 6/iZ 


,a(n 5 / 3 - 1), 


5n 2 


= 3/icr[3n + n 5 / 3 (2n- 


5) + Z|(3n 5 / 3 -5n + 2). , 



(39) 
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Figure 5: Bi-polar electrostatic pulse in real-time units. The solid (red) curves are for highly Lorentzian electrons (k ~ 3/2) 
and the dashed (blue) curves are for Maxwellian electrons. 



The quantities Sri, dni^ — > as n — > 1 and it is easy to see that the function V n ^,\ satisfies all the requirements of 
Sagdeev potential. However, its interpretation is valid only in the limit n — > 1. 

To complete the analysis, we now take the small amplitude limit of the V n -ti in Ea.(|25p by expanding around n = 1. 
The resultant equation, to the third order, is, 



ltd - 2 



■> \ di/ 71 ) +XiSn 2 + X 2Sn 3 = 0, (40) 

which we recognize as the Korteweg de-Vries (K-dV) equation in the transformed space-time co-ordinate 77, describing 
soliton solutions. The factors \i an d X2 are given by 



h v 2 M - M 2 



\, = — 2 Nl) 

1 



X2 = [3/i/ 2 ^ 2 (7a -15v 2 M ) + 9l 2 J. 



2 ( 2 
2 



•4A 2 {9< + 3Z 2 V - v 2 M a(l + 11©}] (42) 



and 



h + he, 



h = /i«m _ /s- 



(43) 
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In the above equations, we note that the factors /i 2,3 > 0. The amplitude A and width A of this K-dV soliton is 
given by, 

A = -^i, (44) 
X2 



a . f A, (45) 

so that for existence of a soliton, xi < 0, which translates to the condition, 

Wf* + <r<v M < V/« + o", (46) 

where the factor 

f = £ = A 2 "- 3 ) (47) 

/K /1 (2k -1) {1-^(1 -/?)}' 1 ' 

represents the effect of two-temperature electrons. Under various limiting cases, we can reduce this condition to simple 
forms, 





1 a < V M < \ z 7z &r + cr, k Y 00 (48) 



1 - v h {l -p) ' V 1 - 14,(1 - P) 



- 3/2) / /3(k-3/2) 

' <y < v M < \ t- 1 — n ' + <r, k^3/2 (49) 



*y i-v h (i-0) \i-v h {i-0) 

which further reduce to (|32l) and (j3"3")) if we consider only single temperature for electrons (v^ — > 1). As the factor 
/« lies between zero and unity, the effect of two different temperatures is not very significant expect when a is small. 
We however note that the parameters /3 and i>h related to two different temperatures can dictate the exact parameter 
regime, where solitons may exist. 



V. CONCLUSIONS 



In this paper, we have discussed the physical situation where electrostatic solitary structures can form in a warm 
plasma immersed in a constant background magnetic field and in the presence of a non-thermal electron fluid. The 
results are viewed in reference to the bi-polar ESWs, typically observed in the auroral regions. We have shown that in 
presence of a super-thermal electron population, large amplitude solitary structure can form in the parallel direction, 
which otherwise is not possible for Maxwellian electrons. We have also emphasised that with a super-thermal electron 
fluid, solitons can form for lower Mach numbers. We have presented our results in the real-time units of the bi-polar 
ESWs, which are basically manifestations of solitary waves, and our results are in good agreement to the available 
experimental data in the auroral region. 

We have further considered the presence of two components of electron populations, both of which can be super- 
thermal, in view of the observations by the certain space-borne experiments viz. Cluster spacecrafts!. However we 
find only marginal modifications to the ESWs by including the two-temperature electrons, and conclude that they 
may have only minor role in dictating the exact parameter details of the ESWs. 
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